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We have studied the ground state of the two-dimensional Hubbard model by using the recently 
proposed adaptive sampling quantum monte carlo (ASQMC) method. We have paid attention to 
the model whose non-interacting band dispersion is almost flat near (7r,0). To minimize the effect 
of the finite size gap overlying on the Fermi level, we have tuned both filling and band structure. 
We found enhancement of the d-wave correlation function at large distance, the spin gap and the 
momentum distribution function consistent with the d-wave gap. We also found the coexistence 
of both the commensurate and incommensurate peaks in S{q), which does not contradict a recent 
experimental finding that both the resonance peak and the incommensurate peaks reside in the same 
doping level of YBCO and BSCCO. 
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It is very difficult to calculate ground state proper- 
ties of large two-dimensional (2D) Hubbard clusters away 
from half filling when C//t is not small because of the no- 
torious negative sign problem for quantum monte carlo 
(QMC) methods. We have recently proposed the adap- 
tive sampling quantum monte carlo (ASQMC) method 
to cope with the problem. The start of the project to 
develop the ASQMC method was influenced by the publi- 
cation of the constrained path (CP) monte carlo method 
by Zhang, Carlson and Gubernatis (ZCG). [|j However 
we found it is possible to reduce the negative sign ratio 
without any constraint if we adopt a new update scheme. 
This may be possible because the update scheme adopted 
by ZCG itself seems to suppress the negative sign ratio, 
which should be the reason why their method is superior 
in accuracy than the more elaborate positive projection 
method proposed by Fahy and Hammann 8 years ago. |^] 
We developed our update scheme within the frame work 
of the " exact update method" of the auxiliary field quan- 
tum monte carlo method. Our scheme is not only easier 
to understand for those who have been working on Hub- 
bard and related models but also it is superior in accuracy 
than the ZCG's scheme because it is free from the inac- 
curate mixed estimator method for measurements and 
its variants which are indispensable to the CP method. 
This should be one of the main reasons why our ASQMC 
method gives more accurate results than the CP method. 
Actually, we observed the ASQMC gives better value for 
the ground state energy of 14 electron 4x4 cluster when 
C//i = 12. Even without constraint, the negative sign 
ratio is reduced largely in the all cases studied, which 
makes measurements stable. < Sign > decrease to 
very slowly as we increase the projecting imaginary time 
T. (Fixing At and Tc constant.) Typically, r = 20 to 
30 is feasible. We have made the serial correlation test 
in many cases and found that samples obtained by the 
ASQMC are statistically independent and the variance 
calculated by the standard formula gives a good esti- 



mate for the true variance. As an example we show the 
autocorrelation function for the d-wave superconducting 
correlation function at (f = in the case of parameter set 
(i) described below in Fig. |^. 

We used the ASQMC method to study the ground 
state properties of the 2D Hubbard model. The parame- 
ter region we paid our attention was such that the non- 
interacting band dispersion was almost flat near the (tt, 0) 
point. We tuned electron filling and band dispersion such 
that the energy difference between the highest occupied 
level (HOL) and the lowest unoccupied level (LUL) was 
close to zero. The Hubbard model we studied may be- 
long to different class than that frequently studied on the 
bipartite lattice. We have studied various parameter re- 
gions of 6 X 6, 8 X 8, 10 X 10 and 12 x 12 clusters, g] 
Here we discuss following three cases; (i) 6 x 6 lattice 
with t' = -0.1667, and t" = 0.2. U = A. The number 
of electron is 28. (ii) 10 x 10 lattice with t' ~ —0.2, and 
t" = 0.05. J7 = 4. The number of electron is 84. (iii) 
12 X 12 lattice with t' = -0.1676, and t" = 0.2. [/ = 4. 
The number of electron is 92. To make sure that our 
result is convergent to the r — > oo limit, we show the r 
dependence of the ground state energy of (i) in Fig. ||. 
The calculation was stable at least up to r = 20, but the 
energy converged at r = 3. We observed similar fast con- 
vergence in the case of (ii). Hence, we proceeded calcula- 
tions with r = 3 to save CPU time. The fast convergence 
suggests that energy gaps among the states with the same 
symmetry are relatively large, while the gap between the 
ground and low lying excited states in the finite size sys- 
tem can be small, if these state have different symmetries. 
This is because we are using the projector method. We 
found the enhancement of the d-wave superconducting 
correlation at large distance in the ground state when we 
used the parameter sets (i) and (ii). We also found the 
momentum distribution function can be well described 
by using the d-wave BCS mean field theory. We found 
exponential decay of the distance dependence of the spin 
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correlation function, which is suggestive of opening the 
spin gap. [Q All these are consistent with the d-wave 
superconducting ground state. Even more surprisingly, 
the commensurate peak coexists with the incommensu- 
rate peaks in S{q) = J dujS{q,uj) of the superconducting 
state. This is characteristic to the parameter region in- 
volving the three parameter sets. The coexistence was 
most clearly seen in a over doped region (iii) which was 
shown in Fig. |^. (The peak interval between the in- 
commensurate peaks is largest in (iii).) Experiments on 
YBCO and BSCCO seem to suggest that both the reso- 
nance peak and the incommensurate peaks exist at dif- 
ferent UJ but in the same filling. Q The suggestion from 
experiments does not contradict our result on S{q). In 
the other parameter regions, we did not so far observe 
similar S{q). The numerical results so far obtained by 
using the ASQMC method suggest that the band disper- 
sion close to (tt, 0) seems to be very crucial to the ground 
state of the 2D Hubbard model. 
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6x6, p=28/36, U/t=4, 
=-0.1667, t"=0.2 
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FIG. 2. The projecting imaginary time r dependence of the 
ground state energy of the 6x6 cluster. 
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